Thermoelectric transport properties of CaMg 2 Bi 2 , EuMg 2 Bi 2 , and YbMg 2 Bi 2 
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The thermoelectric transport properties of CaMg2Bi2, EuMg2Bi2, and YbMg2Bi2 were charac- 
terized between 2 and 650 K. As synthesized, the polycrystalline samples are found to have lower 
p-type carrier concentrations than single-crystalline samples of the same empirical formula. These 
low carrier concentration samples possess the highest mobilities yet reported for materials with 
the CaAl2Si2 structure type, with a mobility of ~740cm 2 /V/s observed in EuMg2Bi2 at 50 K. 
Despite decreases in the Seebeck coefficient (a) and electrical resistivity (p) with increasing tem- 
perature, the power factor (a 2 / p) increases for all temperatures examined. This behavior suggests 
a strong asymmetry in the conduction of electrons and holes. The highest figure of merit [zT) is 
observed in YbMg2Bi2, with zT approaching 0.4 at 600 K for two samples with carrier densities of 
approximately 2 x 10 18 cm -3 and 8 x 10 18 cm -3 at room temperature. Refinements of neutron powder 
diffraction data yield similar behavior for the structures of CaMg2Bi2 and YbMg2Bi2, with smooth 
lattice expansion and relative expansion in c being ~35% larger than relative expansion in a at 
973 K. First principles calculations reveal an increasing band gap as Bi is replaced by Sb then As, 
and subsequent Boltzmann transport calculations predict an increase in a for a given n associated 
with an increased effective mass as the gap opens. The magnitude and temperature dependence of 
a suggests higher zT is likely to be achieved at larger carrier concentrations, roughly an order of 
magnitude higher than those in the current polycrystalline samples, which is also expected from the 
detailed calculations. 



I. INTRODUCTION 

Ternary antimonides with the CaAl2Si2 structure- type 
(AA 2 Sb2; A=Cd,Zn; A=Ca,Yb,Eu) have shown promis- 
ing thermoelectric performance in the 500 to 800 K tem- 
perature range. Thermoelectric efficiency is generally 
discussed in terms of the materials-level figure of merit, 
zT = (a 2 T)/(pn), where a is the Seebeck coefficient, p 
the electrical resistivity, and k the thermal conductivity. 
The highest figure of merit reported within this family 
of materials is zT = 1.2 at 700 K in YbCdi. 6 Zn .4Sb 2 F 1 
while zT values between 0.5 and 1.1 are typicalP^ In 
general, ternary antimonides that are understood using 
Zintl chemistry have received significant attention within 
the thermoelectric community due to their tunability and 
structural complexityP^ However, the analogous bis- 
muthides have received much less attention. 

This work is motivated by electrical resistivity and 
Hall effect measurements on ^4Mg2f3i2 single crystals 
that revealed properties similar to those in ^4Zn2Sb2.^ 
While nominally charge ba lanced and predicted to be 
narrow gap semiconductors} ^ 19 * 20 ! these materials pos- 
sess p-type conduction with carrier densities generally 
on the order of 10 19 -10 20 cm~ 3 , though semiconducting 
behavior was reported for CaCd2Sb2p3 Interestingly, the 



rare-earth containing compounds are reported to have 
higher carrier mobi lity than the alkaline-earth contain- 
ing compounds J2U2I leading to larger zT due to a reduced 
p for a given a. This trend in mobility was observed 
in single crystalline samples of the compounds studied 
here. Transport measurements in the afe-plane revealed 
the single crystals to be heavily defective, and the resid- 
ual resistivity ratios p(300K)/p(10K) were less than 1.6 
for all materials, with the lowest observed for CaMg 2 Bi 2 
(p(300K)/p(10K) = 1.2). These results suggest that de- 
fect scattering significantly affects carrier transport in 
these materials, especially in the alkaline earth com- 
pounds. 

While large single crystals of ylMg 2 Bi 2 were grown 
from nominal melt compositions of AMg^Big, the crys- 
tals were generally irregularly shaped and were not easily 
processed into the geometries required for thermoelectric 
transport measurements at high temperatures. There- 
fore, polycrystalline samples were produced to evaluate 
the thermoelectric efficiency in these compounds. How- 
ever, the as-synthesized polycrystalline samples are found 
to contain different carrier densities than the single crys- 
talline samples. This discrepancy highlights the lack of 
knowledge regarding defects in this family of compounds. 
It is particularly interesting that the polycrystalline sam- 
ples have a lower carrier concentration, suggesting these 
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polycrystalline materials contain fewer acceptor-like de- 
fects, or more electron donating impurities or compen- 
sating defects than the single crystals. 



II. METHODS 
A. Synthesis 

Polycrystalline samples of CaMg2Bi2, EuMg2Bi2, and 
YbMg 2 Bi 2 were prepared via melting, grinding, and an- 
nealing. In a He glove box, the elements were placed 
in a tantalum tube, which was then welded shut under 
an Ar atmosphere; the starting materials were 99.999% 
Bi, 99.95% Mg, 99.9% dendritic Ca from Alfa, and high- 
purity Yb and Eu from Ames Laboratory. The tantalum 
tube was subsequently sealed within an evacuated quartz 
ampoule. Two different melt procedures were employed: 
(1) Melting at 1000°C for 84 h. (2) Melting at 1250°C 
for 30-60 min. Samples produced via (1) utilized stoi- 
chiometric melts, while those synthesized using (2) em- 
ployed 1% excess of Ca, Yb, or Eu. Following melting, 
samples were ground in a He glove box, cold pressed, 
and annealed under vacuum at 600° C for at least 12 h; 
this annealing step took place for all samples made via 
(1) and for the YbMg 2 Bi 2 sample produced via (2). No 
trend between the observed carrier concentration and the 
synthesis procedure was observed. In the following text, 
the synthesis route of a sample is distinguished by num- 
bers corresponding to these procedures, and in the fig- 
ures samples produced via (1) are represented by open 
markers and samples prepared via (2) are represented by 
closed markers. For neutron diffraction measurements, a 
procedure similar to (2) was utilized followed by anneal- 
ing at 600°C for lOh. For procedure (1), total masses 
near 15 g were employed, while 5-7 g was used for (2) and 
for neutron diffraction samples. 

Single crystals were grown from a melt with nominal 
composition AMg4Bi 2 ; here, these samples are denoted 
by the synthesis procedure number 3p3 The excess Mg- 
Bi flux was removed via centrifugation at 650° C, and the 
growth is believed to occur at temperatures lower than 
750°C during which time a cooling rate of « 4.2°C/min 
was employed. Further details can be found in Ref. Qj|] 
Large crystals were obtained with this method, and the 
crystals generally contain one flat face orthogonal to the 
c-axis, and the side opposite this face is typically irregu- 
lar. 

To measure the transport properties, polycrystalline 
powders were hot pressed under dynamic vacuum with 
an applied pressure of ~ 10,000 psi. The graphite fur- 
nace was heated to approximately 650 °C over an hour, at 
which point the temperature was held constant for 1-2 h. 
The samples were hand ground in a He glove box, and 
then loaded into a graphite die using grafoil to separate 
the sample from the graphite die, and remained under 
a He atmosphere until they could be quickly transferred 
to the hot-press. The density of all samples was at least 



97% of the x-ray density. After cutting, samples were 
kept under dynamic vacuum or in a helium glove box 
between measurements. 



B. Crystal structures and phase purity 

Crystal structures were characterized via x-ray and 
neutron powder diffraction. Results obtained are con- 
sistent with those reported in Reference \T5\ for ground 
single crystals. Powder x-ray diffraction (pXRD) data 
were collected at ambient conditions on a PANalytical 
X'Pert Pro MPD using a Cu K a l monochromator. The 
main impurity phase is found to be elemental Bi, with 
volume fractions estimated to be less than approximately 
5% in all samples. 

Neutron powder diffraction measurements were per- 
formed at a series of temperatures ranging from 10 K 
to 973 K, using the POWGEN time-of-flight diffractome- 
ter at the Spallation Neutron Source, Oak Ridge Na- 
tional Laboratory. Samples of CaMg 2 Bi 2 and YbMg 2 Bi 2 
were ground into fine powders in a He glove box, 
and then loaded inside thin-walled vanadium containers. 
EuMg 2 Bi 2 was not studied due to absorption. Measure- 
ments were first performed above 300K, with unsealed 
containers in vacuum inside a radiative furnace. Trace 
oxygen amounts led to a small amount of MgO forma- 
tion. Subsequent measurements below 300K were per- 
formed with the sealed sample container filled with a 
partial atmosphere of helium, and mounted on a closed- 
cycle refrigerator in vacuum. In the configuration used, 
the measurements covered <i-spacings from ~0.3 to 3.5 A. 
A vanadium standard was measured to correct the effi- 
ciency of the raw detector counts.^ Rietveld refinements 
were carried out usin g the GSAS software package and 
the EXPGUI interface 22 23 with (i-spacings from ^0.5 to 
3.0 A. 

Sample homogeneity was investigated using a Hitachi 
TM-3000 tabletop microscope equipped with a Brukcr 
Quantax 70 EDS system. This revealed inclusions of 
MgO in all samples, which were not observed via powder 
x-ray diffraction except for a very small peak in the scan 
of sample Ca2 (Table [I]). MgO inclusions with dimensions 
up to 1=3 50 /im were observed, while the Bi inclusions ob- 
served were generally 1-5 /xm. The extreme broadness of 
the Bi peak in x-ray diffraction suggests the Bi inclusions 
are small and/or strained, and perhaps many are smaller 
than the resolution of this SEM. Oxidation of the surface 
was evident if the samples were not polished immediately 
prior to examination. 



C. Transport measurements 

Thermoelectric transport measurements below room 
temperature were performed in a Quantum Design Phys- 
ical Property Measurement System (PPMS) using the 
Thermal Transport Option (TTO). These measurements 
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utilized gold coated copper leads attached to the sam- 
ples via silver epoxy (H20E Epo-Tek), which requires a 
^30min cure at 100°C. Electrical resistivity and Hall ef- 
fect measurements were performed using a standard four- 
point configuration in the PPMS; samples were thinned 
and electrical contacts were made with 0.025mm plat- 
inum wire and silver paste (DuPont 4929N). Above 
room temperature, the Seebeck coefficient and electri- 
cal resistivity were measured in an ULVAC ZEM-3 M8; 
these samples were cut to rectangular bars between ap- 
proximately 7.5 and 10 mm long and 1.5-2 mm on a 
side. Thermal diffusivity Dt measurements were con- 
ducted on ~10mm diameter discs between 1.5 and 2 mm 
thick, coated with graphite, in an Anter FL-5000 with 
a graphite furnace; data analysis followed ASTM 1461 
for flash diffusivity. The thermal conductivity k (above 
room temperature) was calculated via k = Dt X Cp x d 
where d is the density. The specific heat capacity C p 
was approximated using the Dulong Petit limit C v = 3R, 
calculated from the composition. Based on the equation 
C p = C v + a T T/(d(3T), where ar is the coefficient of 
thermal expansion and Pt is the isothermal compress- 
ibility, in the 300-600 K temperature range, the thermal 
expansion term is small and will cause a slight underes- 
timation of C p . Samples were lightly polished prior to 
measurements . 

The large single crystals utilized to obtain thermoelec- 
tric transport data at low temperatures were not pro- 
cessed into standard geometries, and could not be mea- 
sured at high temperatures. As such, the mean values 
of p for single crystals (see Table |l| were utilized to ob- 
tain the transport dimensions and calculate the electrical 
resistivity and thermal conductivity. The small distribu- 
tion of properties observed in smaller, regularly shaped 
single crystals allows this method to be utilized without 
impacting the trends observed in p and k, though the 
absolute values may be off by up to ~15%. 



D. First Principles Calculations 

Boltzmann transport calculations of the tempera- 
ture and doping dependent Seebeck coefficients were 
performed for CaMg2Bi2, as well as CaMg2Sb2 and 
CaMg 2 As2 for comparison. These calculations are 
based on the electronic structure obtained with the 
all-electron general potential linearized augmented 
planewave (LAPW) method, 1 ^ using the WIEN2k 
codei^SJ Very dense samplings of the Brillouin zone were 
used with well converged LAPW basis sets with local 
orbitals to relax linearization errors and include semi- 
core states.^The zone samplings were done with uniform 
32x32x16 grids, while the cut-offs for the basis set i^ max , 
were obtained with the criterion RK max =9, where R is 
the smallest LAPW sphere radius. The transport cal- 
culations were done within the constant scattering time 
approximation (CSTA), which assumes that the relax- 
ation time is independent of energy, using the BoltzTraP 



codep^as in our previous work on Pb chalcogenides 12MH 
This method gives predictions of a(T) as a function of 
doping with no adjustable parameters by assuming a 
rigid band model and shifting the chemical potential to 
simulate doping. 

The approach for calculating the electronic structure 
is the same as that used in a prior study of halide scintil- 
lators with the TB-mB J functional) 3 -^ and the electronic 
structure of CaMg 2 Bi 2 obtained with this method was 
already reported™ A standard generalized gradient ap- 
proximation (GGA) functional is utilized to relax the in- 
ternal coordinates in the unit cell by total energy mini- 
mization. This was done using the experimental lattice 
parameters from literature^ and the GGA of Perdew 
and co-workers (PBE) P^l The resulting structural param- 
eters were then used to calculate the electronic struc- 
tures using the recently developed modified Becke John- 
son functional of Tran and Blaha (TB-mB J). 33 While 
standard GGA functionals are designed to reproduce to- 
tal energies and structures, but underestimate the band 
gaps of semiconductors,^ the TB-mBJ functional can- 
not be used for total energy calculations but instead 
yields very mu ch im proved band gaps in a wide vari- 
ety of materials! 33 * 35 ^ These calculations were performed 
relativistically, including spin orbit coupling for all ele- 
ments. The LAPW sphere radii were 2.5 bohr for Ca 
and As, 2.6 bohr for Sb, and 2.8 bohr for Bi. The Mg 
LAPW sphere radius was chosen as 2.5 bohr, except for 
CaMg2As2 where a radius of 2.4 bohr was used. 



III. RESULTS AND DISCUSSION 

The polycrystalline samples of AMg2Bi2 produced via 
melting and hot pressing have carrier concentrations be- 
tween 6xl0 18 and ~lxl0 19 holes/cm 3 . While p-type con- 
duction is also observed in the single crystalline samples, 
the free carrier concentration is considerably lower in the 
polycrystalline samples. These results are summarized 
in Table |TJ where standard deviations represent the dis- 
tribution of properties observed for single crystals from 
measurements on multiple crystals as discussed in Refer- 
ence dni Here, the carrier concentration is characterized 
by the Hall carrier concentration, nu = ^/Rhs where 
Rh is the Hall coefficient and e the elemental charge. 

The source for p-type conduction is not known. In 
single-crystalline AMg2Bi2 and polycrystalline AZn2Sb2 
the rare-earth derived compounds have higher carrier mo- 
bilities, and the A=Yb compounds have higher carrier 
densities than the ^4=Eu,Ca compounds. Given this cor- 
relation between A, the carrier density, and mobility, it 
seems reasonable that A vacancies are the primary source 
for the charge carriers. It is interesting that single crys- 
tals grown from an excess of Mg-Bi have higher carrier 
densities, and thus one suspects higher defect levels in 
the single crystals. The formation of A vacancies may 
be promoted by growth in excess Mg-Bi, thus leading to 
larger defect concentrations. Synthesis of polycrystalline 
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TABLE I. Summary of room temperature properties in poly- 
crystalline (P) and single crystalline (S) samples of AMg2Bi2, 
as well as the effective mass calculated using a single parabolic 
band model; uncertainties for single crystal samples come 
from measurements on multiple crystals in Ref. 1191 The sam- 
ple IDs are utilized to identify different samples in the figures 
and the associated number is related to the synthesis proce- 
dure. 
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samples from melts with 1% excess A did not result in 
systematic changes in un- 
it is also possible that Ta incorporates into the poly- 
crystalline samples, which would act as an electron donor 
and thus decrease the hole concentration. While Ta was 
not detected via EDS, the concentration of Ta 5+ nec- 
essary to modify the carrier density is extremely small 
and is well below the resolution of EDS. Finally, growth 
in a sealed vessel may result in higher Mg content, and 
thus Mg vacancies cannot be eliminated as a possible 
source for free carriers. In the related /lZn 2 Sb 2 , the car- 
rier concentrations are generally manipulated via isoelec- 
tronic substitutions, for example in solid solutions based 
on Zn-Cd substitutions or Ca-Yb substitutions. Non- 
isoelectronic doping, such as La for Yb, typically fails 
to change the carrier concentration in these materials. 
Thus, the defects in these materials warrant further in- 
vestigation, and the unintentional manipulation of car- 
rier density in these materials is perhaps one of the most 
important aspects of the current study. 




E 
o 



° 1 



(C) 



0.1 



v ~ 600 

1 w 

> 

cm^ 400 
=f 200 



• m 





Q I i i i I i i i I i i i I 

100 200 300 

T(K) 



A. Low Temperature Transport Properties 

The electrical resistivity p, Hall carrier density tih, and 
Hall mobility (ph=Rh / p) of polycrystalline samples are 
shown in Figure [TJ The temperature dependence of the 
electrical resistivity and Hall mobility are similar in all 
samples, with a suppression of the mobility observed at 
low temperatures leading to an increase in the resistiv- 
ity. This behavior was not observed in single crystals of 
YbMg 2 Bi 2 and EuMg 2 Bi 2 (nor in sample Yb2), though 
sample Ca3 and other CaMg 2 Bi 2 single crystals^ did 
possess qualitatively similar behavior. 

Given the uncertainty in defect identities and concen- 
trations in these materials, one can only speculate on 
the source of decreased pa at low T. The most plausi- 
ble sources are the formation of a mobility edge due to 
disordered defects, the scattering of holes by point de- 



FIG. 1. (a) Electrical resistivity, (b) Hall carrier density, and 
(c) Hall mobility of polycrystalline CaMg2Bi2, EuMg2Bi2, 
and YbMg2Bi2 samples with symbols linked to Table [I] 

fects or ionized impurities, or activation through impu- 
rities (MgO) at the grain boundaries. While activation 
through 'dirty' grain boundaries offers a slightly simpler 
explanation, it does not account for the observation of 
similar behavior in single crystalline CaMg 2 Bi 2 . 

Transport data for single crystals clearly show that de- 
fect scattering strongly influences the mobility in these 
compounds (low resistivity ratios implied by Figure 
|3^) The scattering of carriers by point defects is typ- 
ically modeled with a temperature independent relax- 
ation time, which is generally associated with the residual 
resistance in a material. When point defects are con- 
sidered ionized impurities, however, the scattering typ- 



5 




100 200 300 
T(K) 




T(K) 

FIG. 2. The (a) Seebeck coefficient and (b) thermal 
conductivity in polycrystalline CaMg2Bi2, EuMg2Bi2, and 
YbMg2Bi2 below room temperature, with symbols linked to 
Table [I] The inset in (b) plots k at the low temperature peak 
versus the room temperature carrier concentration, and in- 
cludes data for single crystalline samples. The estimated lat- 
tice thermal conductivity is indistinguishable from the data 
in (b), as the largest estimate for the electronic contribution 
is less than 0.1 W/m/K (for YbM g2 Bi 2 at 150 K). 



ically results in a mobility that increases with increas- 
ing temperature as roughly T 15 P3ln the single crystals, 
defect scattering appears to be stronger in CaMg2Bi2 
than in YbMg 2 Bi 2 and EuMg 2 Bi 2 , and this is consis- 
tent with the trends observed in fiH f° r polycrystalline 
AZn 2 Sb 2 samples at room temperature. In these mate- 
rials, p is increasing with T, which implies that at high 
T phonon scattering limits the carrier relaxation time. 
When acoustic phonon scattering limits the carrier relax- 
ation time, the mobility decreases with increasing tem- 
perature as 7 1-1 ' 5 for non- degenerate carriers and as 
for degenerate carriers. 

The transport properties of CaMg 2 Bi 2 and EuMg 2 Bi 2 
reveal the importance of defect scattering, as a larger mo- 
bility is observed in the low-carrier concentration poly- 
crystalline samples than in the high-carrier concentration 
single-crystalline samples (at low T). To our knowledge, 
these low temperature pu values are the largest mobil- 



ities reported in CaAl 2 Si 2 type materials. The mobility 
in polycrystalline YbMg 2 Bi 2 is suppressed below that of 
the single crystalline material despite a lower defect level 
(Table [I]), which suggests that grain boundary scatter- 
ing may also be important. However, it is important to 
note the irregular temperature dependence of nn for the 
YbMg 2 Bi 2 samples and thus this hh data may be an 
artifact associated with the simple, single band model. 

The thermoelectric transport properties below room 
temperature are shown in Figure [2] for polycrystalline 
samples and Figure [3] for single crystalline samples. The 
lower carrier density in the polycrystalline samples is 
clearly reflected in the increased Seebeck coefficient. The 
Seebeck coefficients decrease with decreasing tempera- 
ture and increasing n#, as expected for simple, doped 
semiconductors. 

The inset in Figure [3^, shows the electrical resistivity of 
EuMg 2 Bi 2 across the antiferromagnetic transition near 
7KP^H While approaching this transition from higher 
temperatures, the electrical resistivity increases slightly 
then decreases sharply below the transition. This be- 
havior is associated with increased magnetic fluctuations 
near the transition leading to increased carrier scattering. 
The Seebeck coefficient also rises sharply and reaches a 
maximum before falling rapidly towards zero. A large 
peak was observed in the heat capacity near this temper- 
ature, and analysis revealed an entropy change consistent 
with the magnetic ordering of Eu 2+ where J = | .^Inter- 
estingly, a similar, yet broad increase in a was observed 
in the single-crystalline CaMg 2 Bi 2 , possibly associated 
with changes in the dominant scattering mechanism. 

The Seebeck coefficient and Hall carrier density can be 
utilized to obtain a density of states effective mass m* 
when a dispersion relationship and a carrier scattering 
model are assumed. Given the decay of the carrier mo- 
bility with increasing temperature above roughly 100 K, 
it is clear the carrier relaxation time is strongly influ- 
enced by acoustic phonon scattering. It is very difficult 
to account for the influence of defect scattering on the re- 
laxation time. Therefore, effective mass values are calcu- 
lated assuming a simple model, with transport in a single 
parabolic band (SPB) that is limited by acoustic phonon 
scattering. This SPB model analyzes the Hall carrier 
density nu = 1/Rhz, which is equivalent ton — rn/Rne 
where the chemical doping level n calculation utilizes the 
Hall factor r# to account for scattering effects that cause 
a difference between n and n#. Further details of this 
SPB model can be found in Reference f3"Tl 

Within the SPB model, the effective masses of holes 
are between 0.2 and 0.6 m e at room temperature, where 
m e is the free electron mass (see Table |l|. For compar- 
ison, the Seebeck and Hall coefficients in AZn 2 Sb 2 are 
well described by a SPB model with m*=0.56m e for 
A=Ca,Sr,Eu,Yb. Some variation between compositions 
exists, as well as between the single crystalline and poly- 
crystalline samples. In addition to possible changes in 
the band structure with composition, the influence of 
defect scattering clearly changes with doping level and 
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FIG. 3. (a) Electrical resistivity, (b) Seebeck coefficient, 
and (c) thermal conductivity in single crystals of CaMg2Bi2, 
EuMg2Bi2, and YbMg2Bi2 below room temperature. In (a) 
and (c) the electrical resistivity and thermal conductivity are 
scaled to match electrical resistivity data in smaller crystals 
(Table [ij to account for the slightly irregular shape of the 
large single crystals. The estimated lattice thermal conduc- 
tivity is within the size of the data markers in (c); the largest 
electronic contribution to k is estimated to be 0.55W/m/K 
for YbMg 2 Bi 2 at 150 K. 



sample composition and this can lead to changes in the 
calculated m* . The validity of a single-carrier model is 
also brought into question for the lowest carrier concen- 
trations due to the appearance of carrier activation at 
only moderately high temperatures, as is the assump- 
tion of a SPB given that electronic structure calculations 



show^ the valence band edge is dominated by two light 
bands (though they possess similar extrema) . A detailed 
theoretical study, such as that presented in Reference [551 
would likely be beneficial for the development of a greater 
understanding of transport in these materials. 

Due to the change in either conduction type or scat- 
tering mechanism at low temperature, the Hall coeffi- 
cient is not analyzed in detail and the Hall carrier density 
n,H=l/ 'Rh I 'e is plotted in Figure [It. At room temper- 
ature, the data suggest r# varies between roughly 1.17 
and 1.15, which is consistent with the low carrier con- 
centrations. A maximum of roughly 1.18 is obtained at 
the non-degenerate limit and thus the 'true' or 'chemi- 
cal' carrier concentrations are likely higher than the Hall 
carrier concentration by about 16%, which does not in- 
fluence the general trends discussed. 

The thermal conductivity data shown for polycrys- 
talline samples (Figure [2k>) and single crystalline samples 
(Figure [3b) reveals that defect scattering also influences 
the thermal transport in these materials. This behav- 
ior is highlighted in the inset of Figure [2]d, which plots 
the value of the low-temperature peak in k versus carrier 
concentration. The decrease in this K max is likely due to 
the point-defect scattering that is associated with hole 
formation. The scattering of phonons by free carriers is 
unlikely in these materials due to low carrier concentra- 
tions and small band masses.^ At higher temperatures, 
phonon-phonon scattering dominates the thermal trans- 
port and k decreases with increasing temperature. The 
influence of defect scattering is observed through a com- 
parison of Figure [2Jd for polycrystalline samples and Fig- 
ure |3j: for single crystalline samples, which reveals lower 
thermal conductivity in the single crystalline specimens 
(in the a6-plane). The single-crystalline data represents 
the high-carrier concentration markers in the inset of Fig- 
ure EJd. 



B. High Temperature Transport Properties 

The high temperature electrical transport properties 
of polycrystalline samples are shown in Figure [4] The 
electrical resistivity first rises with increasing T, due 
to phonon scattering, then reaches a maximum and de- 
creases. The observed maximum is associated with the 
thermal activation of charge carriers, and a correspond- 
ing maximum is observed in the Seebeck coefficients as 
well. 

Assuming charge carriers are activated across the en- 
ergy gap, and not from impurity states, the tempera- 
ture dependence of the resistivity above the maximum 
is expected to be described by p = p Exp[(Eg) / '(2kT)], 
where Eq is the thermal energy gap. While limited 
data exist at these high temperatures, reasonable fits are 
obtained yielding 0.44 eV for Ca2, 0.2 eV for Yb2, and 
0.3 eV for Eul, Eu2, and Ybl. A band gap can also be 
estimated from the maximum in the Seebeck coefficient 
(a m ax at temperature T max ) via E G = 2a max T max eW% 
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For these materials, this relationship yields energy gaps 
of roughly 0.2 eV for the Ca- and Eu-based materials, 
and E G -0.3 cV for YbMg 2 Bi 2 . These estimates of E G 
are somewhat lower than those predicted by first princi- 
ples calculations, which yielded a direct gap of 0.7 eV for 
CaMg 2 Bi 2 P 

As expected from the low temperature data, the poly- 
crystalline samples have large Seebeck coefficients above 
room temperature. The Seebeck coefficients appear rel- 
atively independent of temperature while the electrical 
resistivity decreases with increasing temperature above 
— 450 K for all samples. This is in part due to the scale 
necessary to view all of the data, as the Seebeck coeffi- 
cients generally possess the expected temperature depen- 
dence; a first increases with T then decreases at high T 
consistent with p{T). However, the temperature depen- 
dence of a and p implies strong asymmetry in electron / 
hole conduction. 

The relative rates of decrease in a and p at high tem- 
peratures suggest that hole mobility is much larger than 
electron mobility. This would cause the electronic prop- 
erties to be dominated by hole conduction, and the de- 
creases in a and p would be primarily associated with 
the increase in hole concentration due to the activation 
of electron-hole pairs. The thermoelectric power factor 
{a 2 /p = a 2 a shown in Figure ^) continues to increase 
with increasing temperature despite the appearance of in- 
trinsic conduction, supporting the idea that a large asym- 
metry in electron-hole properties exists. 

In a moderately doped semiconductor with symmet- 
ric band properties, the power factor decreases at high 
temperatures due to the compensation of the Seebeck 
coefficient from minority carriers; the decrease in p is not 
enough to overcome the sharp decrease in a 2 . However, 
deviation from this behavior occurs when the band prop- 
erties are highly asymmetric. In the current scenario, 
electronic structure calculations show that the density of 
states is much larger in the conduction band than in the 
valence bandf^l which is associated with a larger effec- 
tive mass for electrons. The mobility is classically de- 
fined as being proportional to the ratio of the relaxation 
time to the effective mass, and the relaxation time it- 
self is generally a function of the effective mass. When 
acoustic phonon scattering limits the carrier relaxation 
times, the mobility is proportional to m*~ 5 / 2 and the 
conductivity is proportional to to* -1 fRef. [31))), thus an 
increase in the effective mass generally decreases the con- 
ductivity and relative contribution of that carrier to the 
transport properties. Another possibility is that phase 
impurities or impurity states in the polycrystalline sam- 
ples consume the thermally activated electrons. In either 
case, the current data support the claim that hole trans- 
port dominates the electronic transport properties, even 
in the high-temperature intrinsic region. This type of 
behavior is favorable for thermoelectric application. 

The high temperature thermal conductivity data are 
shown in Figure [5^,. A decrease in k with T is ob- 
served, which is expected for these crystalline materials. 




q I i i i i i i i i i i i i i 

300 400 500 600 
T(K) 

FIG. 4. High temperature (a) electrical resistivity, (b) See- 
beck coefficient, and (c) thermoelectric power factor of poly- 
crystalline CaMg2Bi2, EuMg2Bi2, and YbMg2Bi2 samples 
with symbols linked to Table [I] 



However, defect scattering does influence the tempera- 
ture dependence and the decay is not as strong as the 
T _1 expected for phonon-phonon scattering. Additional 
factors influence this temperature dependence, such as 
being near the minimum thermal conductivity. Interest- 
ingly, the thermal conductivity is largest in EuMg 2 Bi 2 , 
and variations in the defect concentration do not trend 
with k at room temperature. Electronic contributions 
are very low due to the large electrical resistivities, par- 
ticularly at room temperature, and thus the differences in 
k between samples are either due to inherent properties 
or variations in sample quality. 

The thermoelectric figure of merit zT rises with in- 
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FIG. 5. (a) Thermal conductivity assuming a heat capacity of 
3RN, (b) thermoelectric figure of merit zT for polycrystalline 
CaMg 2 Bi 2 , EuMg 2 Bi 2 , and YbMg 2 Bi 2 . 



creasing temperature for all samples for all temperatures 
probed, as shown in Figure [HJd. The largest zT is ob- 
tained for YbMg 2 Bi 2 , with both Ybl and Yb2 possessing 
zT -0.35 at 600 K, and zT > 0.4 is obtained by 640 K. 
This value is similar to that reported for isostructural 
Mg 3 Sb 2 , which possesses zT ~ 0.55 at 675 KP The com- 
pound Mg 3 Bi 2 is reported to be metallic] 41 * 42 ! however, 
and thus can be expected to possesses low thermoelectric 
efficiency. 

The sample with higher carrier concentration, Yb2, 
has a larger power factor at lower temperatures, while 
the lower doped Ybl sample has a larger power factor 
at high temperatures due to greater temperature depen- 
dence. This leads to slightly larger zT in sample Yb2 
at lower temperatures, which suggests that higher dop- 
ing levels will result in larger zT and in particular larger- 
average zT. Even the oversimplified SPB model sug- 
gests that sample Yb2 is essentially optimally doped at 
room temperature, and the optimal doping level is ex- 
pected to increase roughly as (m*T) 3 / 2 W^ The optimiza- 
tion of zT clearly needs to be done for performance at 
higher T, as zT is small at 300 K. Therefore, the doping 
levels in these samples are considered low for moderate 
to high-temperature thermoelectric application, and the 
carrier densities found in single-crystalline samples are 
more likely to provide large thermoelectric performance 



at high temperatures. However, the apparent asymme- 
try in electron/hole conduction leads to moderate ther- 
moelectric performance in these samples, and little dif- 
ference is observed in samples with significantly different 
carrier densities. It would be interesting to examine sam- 
ples of the analogous AZn 2 Sb 2 materials with low doping 
levels to see if these trends remain. 



C. First Principles Calculations 

First principles calculations were employed to ex- 
plore the influence of pnictogcn (Pn) on the electronic 
structure and the corresponding transport properties of 
CaMg 2 Pn 2 . This is motivated, in part, by large thermo- 
electric performance in the AZn 2 Sb 2 and ACd 2 Sb 2 com- 
pounds. CaMg 2 Bi 2 , CaMg 2 Sb 2 , and CaMg 2 As 2 are cal- 
culated to be semiconductors, with the calculated band 
gaps decreasing with increasing atomic mass of the pnic- 
tide. The calculated band gaps are given in Table [TT] 
along with the structural parameters that were used in 
the calculation. The gaps are of charge transfer character 
but there is some covalency between Mg and pnictogen 
states. The valence bands have mainly pnictogen p char- 
acter, but with an admixture of Mg, while the conduction 
bands have Mg character, mixed with Pn. These features 
can be observed in the density of states for CaMg 2 Bi 2 , 
which have been published in Reference Q13 These results 
show that the density of states in the conduction band 
is much larger than that in the valence band, as noted 
above, and help to explain the increase in a 2 a with in- 
creasing T at high T. 

The calculated Seebeck coefficients are shown as a 
function of p-type doping level in Fig. [6] at 300 K and 
800 K. The quantity shown is the linear average of the 
c-axis and the two equal in-plane a-axis values. This is 
an approximation to the ceramic average that is made 
because the anisotropy of the electrical and thermal con- 
ductivity are unknown. 

As observed in Figure |6j the As and Sb compounds are 
predicted to have higher values of a for a given doping 
level, even when the doping level is given in units of cm' 3 . 
Since the unit cell volume decreases as one goes from Bi 
to Sb to As, this trend in a(T) is even stronger if one 
considers the number of free carriers per unit cell, which 
is perhaps the more chemical approach. The increase in 
a for a given n is related to an increase in the effective 
mass. This increase in band mass will reduce the carrier 
mobility, though the net effect on the power factor cannot 
be deduced. 

The reduced band gap in CaMg 2 Bi 2 relative to 
CaMg 2 Sb 2 and CaMg 2 As 2 results in a compensation of 
the Seebeck coefficient at high temperatures and low dop- 
ing levels. This is observed as the peak in a near ~3 x 10 18 
holes per cm 3 for T=800K. The influence of bipolar con- 
duction is not predicted until higher temperatures for 
CaMg 2 Sb 2 and CaMg 2 As 2 . Clearly the larger band gaps 
in the As and Sb-based compounds would facilitate op- 
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FIG. 6. Calculated Seebeck coefficient a for CaMg 2 Pn 2 , 
Pn=As,Sb,Bi, for T=300 K and 800 K. The quantity shown is 
the directional average of the Seebeck coefficient as described 
in the text. 



TABLE II. Structural parameters used in the present cal- 
culations. The lattice parameters a and c are experimental 
data from the literature, while the internal coordinates are 
determined by energy minimization. The band gaps, E g are 
obtained with the TB-mBJ functional. 



compound 


o(A) 


c(A) 




ZPn 


E g (eV) 


CaMg 2 As 2 


4.34 


7.13 


0.6333 


0.2481 


1.77 


CaMg 2 Sb 2 


4.66 


7.58 


0.6303 


0.2434 


1.10 


CaMg 2 Bi 2 


4.73 


7.68 


0.6271 


0.2382 


0.70 



eration at higher temperatures. 

Consistent with the above discussion, the magnitudes 
of the Seebeck coefficients suggest that the optimum dop- 
ing level near room temperature is likely to be approx- 
imately lxl0 19 cm -3 , while operation at higher temper- 
atures will require higher doping levels. In the analo- 
gous antimonidcs, high zT is obtained for Seebeck coef- 
ficients between approximately 180 and 250 /xV/K! 1 | 3 | S | 7 1 
In CaMg2Bi2, the calculation results qualitatively sug- 
gest optimal doping levels near 5xl0 19 cm~ 3 for opera- 
tion at 800 K, while the optimum doping is likely to be 
lxl0 20 cm -3 or greater in the As and Sb compounds due 
to the increased band mass. 

The changes in electronic structures imply that the As 
and Sb-based compounds may be interesting thermoelec- 
tric materials. Also, these trends shed light on the way in 
which solid solutions may be used to tune the band prop- 
erties. A detailed investigation into the effects of band 
mass on the power factor is warranted, as is an investi- 
gation of the influence the pnictogen has on the thermal 
conductivity. 



D. Thermal Expansion 

Neutron powder diffraction was undertaken to exam- 
ine the stability and expansion of the lattice. The source 
for free carriers remains unknown, with cation vacan- 
cies being one of the most plausible explanations. In 
YbMg2Bi2, 1 cation vacancy per 1500 unit cells would 
produce lxl0 19 hoies/cm 3 , and thus it is not possible to 
detect these defects via neutron powder diffraction. Sim- 
ilarly, previous single crystalline x-ray diffraction could 
not isolate defects.^ 

Figure [7^, presents the temperature dependence of the 
lattice parameters obtained from neutron powder diffrac- 
tion. The expansion of a and c is similar in both 
CaMg2Bi2 and YbMg2Bi2, as are the values and rela- 
tive changes of the displacement parameters (Fig. [7Jd). 
The refinements possessed R wp of 0.0286 and 0.0282, and 
X 2 of 5.86 and 5.18 for CaMg 2 Bi 2 and YbMg 2 Bi 2 , respec- 
tively, at 300 K. The structural parameters obtained are 
similar to those reported in Ref. Ql|] for ground single- 
crystals. 

Thermal expansion is slightly larger along c for both 
compounds, which is consistent with the first-level ap- 
proximation of an anionic (Mg 2 Bi 2 ) 2 ~ layer being sup- 
ported by charge transfer from the cationic layer of 
A 2+ .^l Figure [7] also reveals these compounds to be stable 
to nearly 1000 K, which is consistent with previous dif- 
ferential thermal analysis to 1073 K. Minor oxidation oc- 
curred during the measurements as observed via the for- 
mation of MgO, the concentration of which increased as 
the measurement temperature increased reaching 2-4% at 
973 K. MgO was also observed in the polycrystalline sam- 
ples after hot-pressing, and the surfaces of hot-pressed 
samples oxidized during the high-temperature measure- 
ments. Thus, while these materials appear thermody- 
namically stable to high temperatures, they are readily 
oxidized and care must be taken to avoid catastrophic 
oxidation, which occurred during one measurement to 
500° C in the Anter system. YbMg 2 Bi 2 survived this 
measurement in the Anter system and appears less prone 
to oxidation than CaMg2Bi2 and EuMg2Bi2- 

IV. CONCLUSION 

The low carrier density observed in polycrystalline 
CaMg 2 Bi 2 , EuMg 2 Bi 2 , and YbMg 2 Bi 2 results in only 
moderate thermoelectric efficiency by 600 K. Two sam- 
ples of YbMg 2 Bi 2 possessed the most promising perfor- 
mance, with zT ~0.4 obtained near 625 K; these samples 
possessed the highest carrier densities of the polycrys- 
talline materials. The thermoelectric efficiency is likely to 
be larger at higher carrier concentrations, with the dop- 
ing levels in the current polycrystalline samples being an 
order of magnitude or more lower than the expected op- 
timum doping levels. First principles calculations reveal 
a decreasing band gap and effective mass as the pnic- 
togen is changed from As to Sb to Bi, behavior which 
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FIG. 7. Refinement results for neutron powder diffraction 
data showing the (a) relative change in lattice parameters and 
(b) the associated change in isotropic displacement parame- 
ters. At 300 K, the lattice parameters were a— 4. 7326(1)A 
and c=7.6780(2)A for CaMg 2 Bi 2 , and a=4.7321(l)A and 
c=7.6649(l)A for YbMg 2 Bi 2 . 



could possibly be utilized to tune band properties and 
enhance thermoelectric performance. The temperature 
dependences of the electrical resistivities and Seebeck 
coefficients suggest that the minority carriers activated 
at high temperatures have little impact on the electrical 
properties, likely due to large minority carrier effective 
mass. Understanding the role of minority carriers and 
the defects responsible for the extrinsic carriers is crit- 
ical for improving the thermoelectric efficiency of these 
and related compounds. 
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